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We analyze the consequences of violation of Lorentz and CPT invariance in the massless neutrino 
sector by deforming the canonical anti-commutation relations for the fields. We show that, for 
particular choices of the deformation, oscillation between massless neutrino species takes place 
when only Lorentz invariance is violated. On the other hand, if both Lorentz and CPT invariances 
are violated, we show that there is no oscillation between massless neutrino species. Comparing 
with the existing experimental data on neutrino oscillations, we obtain bounds on the parameter 
for Lorentz invariance violation. 



I. INTRODUCTION 



There has been an increased interest in the possibil- 
ity that Lorentz and CPT symmetries may be violated 
at very high energies. For example, recent developments 
in quantum gravity suggest that Lorentz invariance may 
not be an exact symmetry at high energies [l[ and CPT 
invariance has also been questioned within such contexts 
0. Spontaneous violation of CPT and Lorentz symme- 
tries can arise in string theories and the violation of 
Lorentz invariance in non-commutative field theories is 
well known [4(. On the experimental side, the UHE (ul- 
tra high energy) cosmic ray events seen at AGASA [E| 
and presently under study by AUGER @ further sup- 
port the possibility that Lorentz and CPT invariances 
may not hold at such energies. Of course, there already 
exist very stringent bounds on Lorentz and CPT viola- 
tion from laboratory experiments in the Kaon and the 
lepton sectors and any violation of these symmetries has 
to be compatible with these limits. Nonetheless, it is 
possible that even a tiny violation of CPT and Lorentz 
invariance can lead to interesting mechanisms for physi- 
cal phenomena. In a recent paper, for example, we have 
shown [7[ how such a violation can lead to baryogene- 
sis in thermal equilibrium (evading one of the criteria of 
Sakharov). In this note, we analyze the consequences of 
Lorentz and CPT violation in the neutrino sector. We 
would like to emphasize that several papers have already 



dealt with the effects of Lorentz [1, [|| and CPT violation 
in the neutrino sector, particularly in connection with a 
qualitative discussion of neutrino oscillation in this sce- 
nario [10] (in another related context see [HI]). In this 
paper, we carry out a quantitative study of such phe- 
nomena within the context of a simple model and derive 
bounds on such symmetry violating parameters from the 
existing experimental results on neutrino oscillation. 

Neutrino oscillation is an interesting phenomenon pro- 
posed about fifty years ago by Pontecorvo (in a differ- 
ent context) which is used to explain the deficit of solar 
and atmospheric neutrinos in fluxes measured on earth 
[H, O EH (for other recent analyses see [IB]). This 
mechanism which is responsible for the resolution of these 
puzzles is closely related to the K°-K° oscillation [l6| ■ In 
its simplest form, the probability for oscillation between 
two species of particles i,j with a mixing angle Oij and 
energy levels Ei , Ej in a time interval t is given by 



p i-*j(t) = sin 2 (2%) sin* 



V 2 



where 



AEij — Ei — Ej. 



(1) 



(2) 



If the oscillation is between two neutrino species z/£, Vj 
with small masses uii,mj respectively, then in the con- 
ventional scenario one expands (this assumes Lorentz in- 
variance and c = 1) 
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so that we have 
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AEij = Ei — Ej 



P+^Z=P+7±, (3) 
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where we have assumed that for neutrinos of small mass, 
Ei ss Ej = E. In this case, the probability for oscillation 
between the two neutrino species in traversing a path 
length L can be written as (see ([T]) and Igj)) 

Pu^iL) = sin 2 (2%) sin 2 I \ 

„ , x , / I.27Am 2 L\ , 
= sin 2 (2%) sin 2 f j , (5) 

where Amfj = mj — m 2 is taken in (cV) 2 , the neutrino 
energy E in MeV and the length of path traversed in 
'm' (meters) (In the last line of the above formula, we 
have restored all the nontrivial constants as well as traded 
the time interval for the path length assuming that the 
neutrino travels almost at the speed of light.) 

It follows from eq. ([5]) that neutrino oscillation does not 
take place in free space if neutrinos are massless or (when 
massive) are degenerate in mass. With three families of 
neutrinos, there can only be two independent combina- 
tions of squared mass differences, say Am 2 2 , Am| 3 which 
are sufficient to find a solution for the solar neutrino as 
well as the atmospheric neutrino puzzles. Within the 
standard model, this can be achieved with the bounds 

m 

Am 2 2 < l(T 4 eV 2 , lCT 3 eV 2 < Amj 3 < 10~ 2 eV 2 . (6) 

Given these, the bound on Am 2 3 = Am 2 2 + Am^ is 
determined. There is no further freedom within a model 
with three families of neutrinos. 

Several experiments by now have looked for neutrino 
oscillations. One such experiment, namely, the LSND 
(Liquid Scintillator Neutrino Detector at Los Alamos) 
|18| has used muon sources from the decay tt + — > /i + + 
v^. The experiment looks for neutrino oscillation in the 
subsequent decay of the muon through /i + — * e + +i/ e +9 ll . 
After a path length of L = 30m, the experiment finds 
the oscillation channel — > v e (with 20MeV < E v < 
58.2MeV) with a probability of 0.26%. The experiment 
also reports the existence of the oscillation — > v e with 
the same probability. Furthermore, the analysis of the 
results of this experiment, following ([5]), leads to a bound 
on the difference of the relevant squared mass difference 
to be 

Am 2 < leV 2 . (7) 

However, the MiniBoone experiement, which was ex- 
pected to verify the results of LSND, has recently re- 
ported their first result [lj| and excludes the mass region 
in 0. As a result, the simple explanation of the LSND 
results, based on the two flavor neutrino oscillation is 
ruled out. 

All of the above discussion has been within the con- 
text of the standard model with a massive neutrino where 
both Lorentz invariance and CPT are assumed to hold. 
On the other hand, if Lorentz invariance or CPT or both 



are violated in the neutrino sector, it has been suggested 
that neutrino oscillation can take place in free space even 
for massless neutrinos (in contrast to eq. ([5]) where 
Lorentz invariance is assumed). This was pointed out 
by Coleman and Glashow [2(| and developed more ex- 
tensively by Kostelecky and collaborators [2ll . |22| . This 
is particularly clear from eq. ([T]) where we see that 
the probability of oscillation really depends on the dif- 
ference in the energy of the two neutrino species and 
if AEij — Ei — Ej 7^ even when the masses vanish, 
the probability of oscillation will be nontrivial. This can 
happen, for example, if the two neutrino species have dif- 
ferent (energy) dispersion relations. This possibility has 
been discussed extensively in the last few years by var- 
ious groups [2l], [23[. In particular, ref [2l| analyzes the 
structure of the most general Lagrangian with violations 
of Lorentz invariance and CPT in an attempt to un- 
derstand the discrepancy between solar and atmospheric 
neutrinos and the LSND anomaly [l8j . 

The goal of this note is to analyze the consequences of 
violations of Lorentz and CPT invariances in the mass- 
less neutrino sector from the point of view of a non- 
commutative field theory where such violations are more 
natural. In such a model, the violation of Lorentz and 
CPT invariances is implemented through a deformation 
of the canonical anti-commutation relations for the neu- 
trino fields. Such a model can be thought of as a subclass 
of the general model proposed in [21| , but since, depend- 
ing on the deformation, we have fewer arbitrary param- 
eters, we naturally have more predictive power. The re- 
sult of our analysis can be summarized as follows. If 
there is violation of only Lorentz invariance, then oscil- 
lations between massless neutrino species can take place 
and comparing with the existing experimental data, we 
can determine bounds on the parameter characterizing 
Lorentz invariance violation. On the other hand, if both 
Lorentz and CPT invariances are violated, there is no 
oscillation between massless neutrino species. 

II. THE MODEL AND THE 
PHENOMENOLOGY OF NEUTRINO 
OSCILLATION 

The model that we will describe below is inspired by 
the quantum theory of non-commutative fields devel- 
oped in [24l ]. The quantum theory of fermionic non- 
commutative fields (neutrinos) is obtained from the stan- 
dard fermionic quantum field theory by deforming the 
anti-commutation relations while retaining the usual 
Hamiltonian. In order to explain in some detail the con- 
struction, let us consider the conventional Lagrangian 
density for two flavors of massless fermions (neutrinos) 
given by 

L = if^d^\ (8) 

where the superscript % = {1,2} runs over the flavor 
quantum number (sum over repeated indices is under- 
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stood). 

The Hamiltonian density has the form 

H = -i U^S • , (9) 

where a — 7°7- With the conventional canonical anti- 
commutation relations for the fermion fields, one would 
obtain the standard relativistic equations for the mass- 
less neutrinos using the Hamiltonian following from ([5]). 
However, the non-commutative theory is obtained by de- 
forming the canonical anti-commutation relations while 
maintaining the form of the Hamiltonian density @ . 

We postulate the deformed equal-time anti- 
commutation relations to have the form (with all 
others vanishing) 

{^(x),^ t (y)} = ^^(x-y), (10) 

where a, (3, ... = 1,2, 3, 4 are spinor indices and A^g is a 
constant matrix. In this paper we consider the following 
two special choices for the structure of the deformation 
matrix A. 

1. A has a nontrivial structure only in the flavor space. 

2. A depends on both the flavor and the spinor indices 
nontrivially through a constant background vector. 

As we will see, the first choice leads to a violation of 
Lorentz invariance whereas both Lorentz and CPT in- 
variances are violated with the second choice. 

1. Deformation depending only on flavor indices 

In this case, A is a 2 x 2 constant matrix with complex 
elements in general, which, for simplicity, can be chosen 
to have the form 

a- = ( i ? ) . (id 

so that the complex parameters a can be thought of as 
the parameters of deformation. Clearly, the deformed 
anti-commutation relations reduce to the conventional 
ones when the parameters of deformation vanish. 

Given the deformed anti-commutation relations (fTO]) 
and the Hamiltonian density , the dynamical equation 
takes the form 

ft = -A ij (a-Wijp) , (12) 

which in momentum space takes the form 

Eft =A ij (a-p ft) . (13) 

In order to determine the energy eigenvalues for this sys- 
tem, let us consider the unitary matrix D 

D i = = -L ( H f). (14) 



It is straightforward to check that D diagonalizes A and 
as a result, the energy spectrum for the fermions follows 
to be (c = 1) 

El = ±(1 + H)|pl, 

El = ±(l-\a\)\p\. (15) 

Here E 1 ' 2 are the energies of the two species of the neu- 
trinos considered. We, therefore, conclude that with this 
choice of the deformation, this system exhibits violation 
of Lorentz invariance, as can be seen from the dispersion 
relations (fT5|) . However, CPT symmetry remains intact 
in this case which can be easily seen as follows. 

As we have emphasized, our deformation can be 
thought of as a subclass of the standard model exten- 
sion (SME) [2j . Indeed, if we restrict to the part of the 
Lagrangian in Q given by 

C = ift^d^ft + iftc^^d^ft , (16) 

with the constant background field cf*"^ diagonal in the 
space-time indices, then we obtain n l = i{8 1 ^ -\-(p°> % i ')ft^ . 
This leads to the canonical anti-commutation relations 

{^(x),^'t(y)} = (^+c 00 «)- 1 5(x-y), 

= A ij 8(x-y), (17) 

which is in agreement with our deformation for the non- 
commutative fields. In this case, since the background 
field is a constant second rank tensor, the extra term 

ft^^d v ft, 

violates Lorentz invariance but not CPT symmetry as 
was also pointed out in reference [2j. 

The energy eigenstates corresponding to the eigenval- 
ues (|15p can now be determined directly through the ap- 
plication of the diagonalizing matrix 




where ip 1 and ip 2 are eigenstates with energy values E 1 
and E 2 respectively. 

The time evolution for the energy eigenstates, is deter- 
mined to be 

ft(t) = e - iE + t+i ^ S ^(0), (19) 
^, 2 (t) = e - tE + t+l P- s yl 2 (0), (20) 

where E^ 2 are the energy eigenvalues determined in (|15|) . 
The diagonal wavefunctions can be seen from (fT4|) and 
(|18p to have the explicit forms 

v-M-^*)- (21) 
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Thus, we see that if a S 5i, then the diagonal wavefunc- 
tions can be thought of as resulting from a rotation by 
7r/4 in the flavor space. This is, in fact, consistent with 
the hypothesis of large mixing angle (LMA) [26| and, 
therefore, for simplicity let us choose a to be real. In 
this case, we can parametrize (f2Tj) as 



(22) 
(23) 



Tp 1 = cosOuip 1 + sm9i2ip 2 , 
ip 2 = - sin6>i2'0 1 + cos^V' 2 



with the mixing angle 812 = 45° . 

Relations (J22) and ([25)) can now be inverted to give 



xp 1 = cos 6'i2'0 1 — sin Oi 2 ip 2 , 



4' 2 



sin 0i2 



■p 1 +COS6>12'0 2 



(24) 



Thus, a neutrino initially in the state ip l would evolve in 
time as 

ijj l (t) = cos6»i2 ^(t) - sin6»i2 %b 2 (t) 

-iE]t 



cms" W12 e + sin 2 0i2 e _lB +* ) ^(0) 



sin20 12 (e~ lB +' 



-iB±t 



^ 2 (0) 



Therefore, at a later time t, the probability of finding the 
state ip 2 in the beam is given by 



isin20 12 (e~^ 4 



= sin 2 (20i 2 ) sin 2 (|a||p1t), 



(25) 



and since we are considering particles with velocities close 
to c, we can replace 

t -> L, 

where L denotes the path length traversed by the neu- 
trino. Let us note here that in our theory, the velocity 
of the neutrino can in principle be different from c, but 
any further correction is suppressed by terms of the order 
0(a 2 ) which is extremely small. Thus, the probability 
for oscillation (l25l) becomes 



(20i 2 ) 



E L) 



(26) 



where we have used the fact that for \a\ <C \p\. 

There are several things to note from the expression 
([26]) . First, the deformation parameter a leads to rota- 
tions in the flavor space and thereby determines the mix- 
ing angles. However, differences from the conventional 
description of massive neutrinos arise because these de- 
formation parameters also determine the nontrivial dis- 
persion relations for the energy eigenvalues and lead to a 
nontrivial energy difference even in the absence of masses. 
Consequently, oscillation takes place even for massless 
neutrinos. The difference from the conventional descrip- 
tion of neutrino oscillation shows up in Q26[) in the fact 
that the energy dependence is linear as opposed to the 
inverse dependence in 



2. Vector- dependent deformations 

In order to define the minimal matrix A which breaks 
both CPT and Lorentz symmetries and which connects 
different flavors -that is A 12 ^ for two flavor indices - 
we need to include a constant background (real) vector 
e]i . With this, the minimal A has the form 



(27) 



If we assume rotational invariance, we can set the space 
components of e^f to zero. The simplest case mixing fla- 
vors would then correspond to choosing ej 2 — e — eg 1 
and the equations of motion in this case would have the 
forms 



■ijj 1 = -a ■ Vifi 1 - ej ■ Vtp 2 , 
?P 2 = -a- V^ 2 - ey- V^ 1 - 



(28) 



These can be diagonalized and lead to the dispersion re- 
lations 



El 



El = ±vT 



\P\, 
\p\- 



(29) 



This is identical to the conventional dispersion relation 
for massless neutrinos, expect for a scale factor, and 
shows in particular that massless neutrinos and anti- 
neutrinos of different species are degenerate in energy. 
As a consequence, oscillations between massless neu- 
trino species cannot take place (see eq. (]2j) if Lorentz 
and CPT symmetries are simultaneously violated in this 
model. 

Equations (|28|) can be also obtained from the Hamil- 
tonian density 



(30) 



with canonical anti-conmutation relations instead of (JTT 
Our deformation can also be understood as a subclass of 
the standard model extension of [2| as follows. Let us 
consider the part of the extended model of the form 



(31) 



with our minimal choice for the constant background vec- 
tor . In this case, the derivation of the canonical anti- 
commutation relation leads to the deformation discussed 
above up to a field rescaling. Since the constant back- 
ground field is a vector, the extra term in the Lagrangian 
violates CPT invariance and the Lorentz violation of this 
term is manifest as well. Thus, the second choice of the 
deformation violates both Lorentz and CPT invariances. 
The CPT violation can also be seen from the analysis 
of the energy eigenstates. The diagonalized fermions in- 
volve a combination of left and right handed fields and, 
as a consequence, do not have well defined CPT trans- 
formation properties. 
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3. Bounds for a (Lorentz invariance violation) 

As we have shown, oscillation between massless neutri- 
nos can take place in our model if only Lorentz invariance 
is violated. Therefore, let us use the existing experimen- 
tal data on neutrino oscillations to derive bounds on the 
parameter a of Lorentz invariance violation (LIV). We 
note that our earlier analysis for two neutrino flavors can 
be extended to incorporate more flavors easily. For ex- 
ample, to accommodate three neutrino flavors, we need 
to generalize the deformation parameter (as well as the 
mixing angles) as 

a—>onj, 612- *0ij, (32) 

where i,j = 1, 2, 3. In this case, the probability for oscil- 
lation between neutrino flavors can be written as 

P Vi ^ Vj {L) = sin 2 (2%) sin 2 (\a t] \EL) . (33) 

The same formula also holds for anti-neutrino oscillations 
since CPT is not violated in this case. 

Therefore, comparing with the conventional analysis 
of oscillation for massive neutrinos given in |5J , we can 
identify (c = 1) 



We note here that there are three deformation parame- 
ters Oiij, without any further constraint unlike the con- 
straint on the difference of the squared masses in the 
conventional scenario. 

Let us next note that from the solar neutrino experi- 
ments, we know that this involves oscillations of the fla- 
vors 1 — ► 2 with 

Am 2 2 < 8 x 10" 5 eV 2 , E ~ lMeV. (35) 

From (|34p , this translates into a deformation parameter 

|ai 2 | < HT 17 . (36) 

The atmospheric neutrino results, on the other hand, in- 
volve an oscillation of the type 2 — > 3 with 

Am 2 ^ < 2 - 6 x l(T 3 eV 2 , E ~ lGeV. (37) 

From (|34p. we see that this would translate into a defor- 
mation parameter 

|« 23 | < 10- 22 . (38) 

Finally, we note that although the simple interpreta- 
tion of the LSND results has been disproved by the Mini- 
BooNe experiment, it is nonetheless interesting to recog- 
nize that here the oscillations involve flavors of the type 
1 — > 2 (both in the neutrino as well as the anti-neutrino 
channels) with 

Am 2 2 < leV 2 , E ~ 50MeV. (39) 



In this case, the analog of (|34|) for the antineutrinos leads 
to 

|a u | < 1(T 16 . (40) 

It is clear now that within this scenario, all the exper- 
imental results can be naturally explained without any 
particular puzzle. We would like to note here that our 
analysis for solar neutrinos are not in contradiction with 
the data from the KamLAND experiment [27| ■ 

In this discussion, we have assumed the neutrinos to 
be completely massless in which case, the conventional 
oscillation does not take place. It is possible that the 
neutrinos have a small mass and that both mechanisms 
do contribute to the phenomena of neutrino oscillation. 
In this case, a careful analysis of the atmospheric neu- 
trino oscillation results can lead to even a more stringent 
bound on the parameter «23 [HI- 

III. CONCLUSION AND OUTLOOK 

In this paper, we have carried out a quantitative anal- 
ysis of the consequences of CPT and Lorentz invariance 
violation in the massless neutrino sector. While it has al- 
ready been suggested that in such a case, neutrino oscil- 
lation can take place even for massless neutrinos, we have 
presented a simple model of a theory of noncommutative 
fcrmions to study this phenomenon quantitatively The 
model contains a minimal number of symmetry violating 
parameters that are introduced as deformation parame- 
ters in the equal-time anti-commutation relations for the 
fermion fields. Real values of these deformation param- 
eters naturally lead to the large mixing angle scenario. 
While the deformation parameters directly lead to mixing 
between different neutrino flavors, they also lead to non- 
standard (energy) dispersion relations (through Lorentz 
and CPT violation), which leads to oscillations between 
massless neutrino species if only Lorentz invariance is vi- 
olated. 

In the case that there is violation of only Lorentz in- 
variance, we have determined bounds on the parameters 
of deformation (parameters characterizing Lorentz invari- 
ance violation) from the existing experimental data on so- 
lar and atmospheric neutrinos as well as from the LSND 
data. The bounds on the deformation parameters within 
this minimal model are obtained to have the values 

|a 23 | < 1(T 22 , |« 12 | < 1(T 17 (41) 

We note that bounds for a were obtained from the LMA 
scenario by using 9 ~ tt/4. For the solar and atmospheric 
neutrinos, this is indeed consistent with the experimental 
determination. 

On the other hand, if Lorentz and CPT invariances are 
violated simultaneously, then oscillation between mass- 
less neutrinos disappears in our model, although oscilla- 
tions can take place for massive neutrinos in the standard 
scenario. This fact could indicate that at very high en- 
ergy where masses can be neglected, neutrino oscillation 



6 



would signal a violation of Lorentz invariance and not of 
CPT symmetry. 
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